ORTHOGONAL POLYNOMIALS 
IN SEVERAL NON-COMMUTING VARIABLES. II 



T. BANKS AND T. CONSTANTINESCU 



Abstract. In this paper we continue to investigate a certain class of Hankel- 
like positive definite kernels using their associated orthogonal polynomials. 
The main result of this paper is about the structure of this kind of kernels. 



1. Introduction 



Positive definite kernels are studied for their manifold applications. In this pa- 
per we consider a special type of kernels K defined on the free semigroup on N 
generators with the property that 



for any words a,a,T, where I {a) denotes the word obtained by writting a in the 
reverse order. These kernels appear in many situations, see for instance [HI and [Jj. 
Our goal is to determine an explicit structure of the positive definite kernels satis- 
fying the above invariance property. Since in the case iV = 1 such kind of kernels 
are precisely the Hankel kernels, it is quite natural to consider associated orthog- 
onal polynomials and to study their properties. Our main result establishes the 
connection between moments and Jacobi coefficients, as a multivariable extension 
of a classical result. We also describe the Jacobi coefficients of the free products of 
orthogonal polynomials. 



We introduce orthogonal polynomials on several hermitian variables and we dis- 
cuss several general results. Especially, we emphasize the usefulness of a matrix 
notation that reduces very much the degree of complexity and makes clear the 
analogy with the classical, one-dimensional case. Let be the unital free semi- 
group on N generators 1, . . . , with lexicographic order -<. In particular, Fj*" is 
the set No of nonncgativc integers. The set of positive integers will be denoted by 
N. The empty word is the identity element of F^ and the length of the word a is 
denoted by \a\. The length of the empty word is 0. There is a natural involution 
on given by . . .ii) = i; . . . ii as well as a natural action of F^ on itself by 
juxtaposition, (tr, r) — > ctt, ct, t G F^. Let Vn be the algebra of polynomials on 
N non-commuting indeterminates Xi,. . .,X]y with complex coefficients. For any 
(T = «i • • • i/ g F^, we define = Xi^ ■ ■ ■ X-i^ . Using this notation, each element 
P E Vn can be uniquely written as 



K{aa,T) = K{aJ{a)T) 



2. Orthogonal polynomials 



(2.1) 




with only finitely many coefficients Ca different from zero. The length of the highest 
a such that Cq. ^ is the degree of P. We also have 

^ = E^^ = E E 

k>0 k>0\a\ = k 

where each Pk belongs to the vector space £^ of homogeneous polynomials of degree 
k > in N variables Xi, . . ., Xjy. The dimension of £^ is N'^. An involution + 
on Vn can be introduced as follows: X'^ = Xk, k — 1,...,N; on monomials, 
(Xa-)^ = Xji^^y, in general, if P has the representation as in H2.1|l then 

(2.2) ^+ = E ' 

and Vn is a unital, associative *-algebra over C. Let (f> he a, strictly positive 
functional on Vn, that is, (/> is a linear unital map on Vn and (t){P^P) > for 
every P G Vn — {0}. The Gelfand-Naimark-Segal construction applied to gives a 
Hilbert space Ti.^ such that {X„}^^f+ is a linearly independent family in 7i^. The 
Gram-Schmidt procedure gives a family {1^90,}^^^+ of polynomials such that 

(2.3) V'q = ^ aa^fjXp, Qa^a > 0; 



(2.4) {<Pa,'Pi3)^ ^ 6a,f3, a,/3eF^, 
where for Pi, P2 G 7^ at. 

The elements (^q, a G F^, will be called the orthonormal polynomials associated 
with (j). We notice that the use of the Gram-Schmidt process depends on the 
order chosen on F^. A different order would give a different family of orthogonal 
polynomials. Due to the natural grading on F^ it is possible to develop a base 
free approach to orthogonal polynomials. In the case of orthogonal polynomials 
on several commuting variables this is presented in ,8j. However, in this paper we 
consider only the lexicographic order on F^. 
The moments of (j) are 

Sa = HX^), a e F+, 

and we define the moment kernel of (j) by the formula K^{a, (3) — s/(q)/3, a, (3 G F^. 
We notice that is a positive definite kernel on F^ and for a,a,T g F^, 

(2.5) K^{aa,T)^K^{<j,I{a)T). 

This property can be viewed as a Hankel type condition. Conversely, it is easily 
seen that if a positive definite kernel K satisfies 1)2. 5|) then there exists a positive 
functional (j) on Vn such that K = Ka,. 



2.1. Three term relations. Let (/> be a unital strictly positive functional on Vn 
and let {v'a}„g][r+ be the orthonormal polynomials associated with 0. As in the 
commutative case (see 0), it is very convenient to use a matrix notation, <!>„ = 
[(^o.]|CT|=n for n > and $_i = 0. With this notation, the analogy with the classical 
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case = 1 will be much more transparent. It turns out that the family {$„}n>o 
satisfies a three-term recursive formula, 

(2.6) Xfe$„ = $„+iA„+i,fe + $„B„,fe + $„_iA; 

for fc = 1, . . . , iV and n > (see and Each matrix Bnj^, n > 0, k — 1, . . . , N, 
is a selfadjoint iV" x N" matrix, while each An^k, n > 0, k — 1,...,N, is an 
N"- X 7V"~^ matrix such that 

An = [ An,l ■ ■ ■ An,N ] 

is an upper triangular invertible matrix for every n > 0, with strictly positive 
elements on the diagonal. The fact that An is upper triangular comes from the 
lexicographic order that we use on F+ . The invertibility of B„ IS a consequence 
of the fact that </> is strictly positive and appears to be a basic translation of this 
information. The diagonal of An is strictly positive since we chose aa^a > 0. A 
family A — {An^k, Bm,k \ n > 0,m > 0, k = 1, . . . ,N} of matrices satisfying all 
these properties will be called admissible. It turns out that there are no other 
restrictions on the matrices An^k, Bn,k as shown by the following Favard type 
result mentioned in j^. A similar result for the monic orthogonal polynomials, 

Pa — If a was recently mentioned in 

Theorem 2.1. Let Lp^ = 'Ylr^a^'^:'^-^'^ ' ^ ^JV' elements in Vn such that 
ifijj — I and flo-^o- > 0. Assume that there exists an adm,issible family A of matrices 
such that for k — 1, . . . , N and n > 0, 

where $„ ~ [^a-]\a\^n f'^''^ n > and <I>_i = 0. Then there exists a unique strictly 
positive functional <f) on Vn such that {<y3(j}o.gF+ family of orthonormal poly- 

nomials associated to cf). 

There is also a family of Jacobi matrices associated with the three-term relation 
in the following way. For P S Vn define 

Since the moment kernel has the Hankel type structure mentioned above in H2.5|l . 
it follows that each (P) is a symmetric operator on the Hilbert space Ti^ with 
dense domain Vn- Moreover, for P,Q £ Vn, 

^^PQ) - *4P)*4g), 

and '^ci,{P)T> C V, hence is an unbounded representation of Vn- Also, 0(P) = 
(^"0(^)1, 1)0 for P e Vn- We distinguish the operators '^k = ^^{Xk), k = 
1, . . . , TV, since *0(I]creF+ c<tXc,) = J2ae¥% c<t*<j. Let {ei, . . . , cn} be the stan- 
dard basis of and define the unitary operator W from P{¥^) onto Hcf, such 
that W{ea) — fcr, o- £ F^. We see that W~^V is the linear space Vq generated by 
Bct, cr e F^, so that we can define 



Jk = W~HkW, A: = l,...,iV, 
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on Vq. Each Jk is a symmetric operator on Vq and by 
closure of) Jk with respect to the orthonormal basis {eo-} 



the matrix of (the 



o-eF 



+ IS 



Jk 



O.k 



A* 



i,fc 

^i.fe Bi,k A. 







2,fc 







2,fe 



We call (Ji, . . . ,Jn) a Jacobi A^- family on Pq- It turns out that the usual admis- 
sibility conditions on A„^k and B„^k insure a joint model of a Jacobi family in the 
following sense. 

Theorem 2.2. Let ( Ji, . . . , J^) be a Jacobi N -family and assume that the corre- 
sponding A is an admissible family of matrices. Then there exists a unique strictly 
positive functional <j) on Vn with associated orthonormal polynomials {</'o-}crgF+ 
such that the map W{e„) = if„, a G F^, extends to a unitary operator from Z^(F^) 
onto and Jk = W-^^^kW for k = I, . . . , N . 

For details about the proof of this result see [5]. 

2.2. Jacobi N -families and combinatorics of lattice paths. The matrices An^k and 
Bn^k contain the whole information about the orthonormal polynomials (or the 
moment kernel K^). Ususally they are called the Jacobi coefficients of and can 
be calculated from the moments. For instance, 



A — \ 1^-^ \ 1®^ 

— Fq,/3J|q| = |/3|=„ [^Q,/3J|q| = |/3|=„ 



-1 ' 



where a^^^ are the coefficients of the orthogonal polynomials and for a matrix A 
we use the notation 



^9 



Ai 



.A. 



I times 

In their turn, the coefficients aa^p, (3 ^ ol, can be calculated from the formula 



\J Da-lDa 



where 



D„-det [K{a\l3')]^,j,^^. 

The formula for Bn.k is somewhat more involved and we do not record it here. 
Instead we consider a different kind of relation between moments and Jacobi coef- 
ficients which appears to be more exphcit. The case = 1 is classical, see [Hj, [TU] . 
The Jacobi iV-family (Ji,...,Jjv) is a convenient tool to deal with this matter. 
Thus, for any a € we have that 



and by Theorem 12. 21 



(/)(X,) = (*4X,)1,1)^ = (*,1,1)^, 



therefore we have 



(2.7) s,^ = (J<,eo,eo), a eWj^. 

Now we introduce some special paths on No x {1, . . . , A^} x Nq. The allowed 
steps are the following: level steps l^ „^ from a point (n, fc, m) to (n + 1, k, m), level 
steps Ip"^ f. from a point (n, k, m) to (n,p, m), for some p e {!,..., N} — {k}, rise 
steps „j from a point (n, fc, m) to {n + 1, A;, m + 1), and fall steps from a 

point (n, k, m) to (n + 1, fc, to — 1) (see Figure 1 for an example). 




Figure 1 . An example of a path ior N = 2 



We introduce a weight on steps by the formula 



t«(step) = < 



where / denotes the identity matrix of appropriate size. If p is made of I steps, 
step 1, . . ., step I, then we define the weight of p by the formula 

u;(p) = w(step I) . . . ■w;(step 1). 

Any word a € — {0} has a unique representation a = i^^ . . . ip'' with li, . . ., 
ip G {1, . . . , N}, ki, . . ., kp > 0, and ii ^ ij+i for Z = 1, . . ., p — 1. We consider the 
set M.a of all paths that start at (0, ip, 0) and end at (|(t|, zi, 0), with the property 
that the first kp steps belong to No x {ip} x No, the next fcp_i steps belong to 
No X {ip-i} X No, and so on, until the last k\ steps which belong to No x {ii} x No- 
These sets are related to the set of Motzkin paths. The Motzkin paths of length n 
are the paths in Nq made of level, fall, and rise steps, starting at (0, 0) and ending 
at (n, 0). Their set is denoted by M.^ and the number of elements of is given 
by the Motzkin number 

h 

It is easily seen that for any a G F+ - {0} there is a bijection between M.^ and 
M.\a\. We can now describe a combinatorial structure of the moments. 

5 



/ if step = 
Bm,k if step - 
Am+i,k if step 



if step = 



Theorem 2.3. Let (j) he a strictly positive Junctional on Vn md let A be the 
admissible family of matrices associated with (j) by (|2.t)|) . Then the moments of (j) 
can be calculated by the formula 

(2.8) E "'(P)' ^eF+-{0}. 

Proof. We consider the following points in Nq x {!,..., N} x Nq: 
F„,, = (0,j, n), n>0,j = l,...,7V; 
QnM,m ^ {n,k,m), m,n> 0,k = 1, . . . ,N. 
For CT e F+ - {0}, CT = i^'i . . . ip" , we claim that 



Jrr - 



•^0.0 ^0,1 
jty ja 
•^1,0 '^1.1 



where the entry ■ gives the sum of (the weights of) the paths in Nq x {1, . . . , iV} x 
No from -P/.i^ to Q\a\.ii,k- The claim is clearly true for \a\ = 1 and then suppose it 
true for any word of length < n. Then consider a word a of length rt + 1. Several 
cases can occur. 

Case 1. k = 0. Let a = i^^ . . .ip'' = iiT. First, assume fci = 1. Due to the fact 
that the level steps of type Ip have weight / and by the induction hypothesis, we 
deduce that the sum of the paths from Pj.ip to Q\a-\,ii,o is 

which is precisely the (0, j) entry of the product 

The case fci > 1 is similar, just by the induction hypothesis we deduce that the 
sum of paths from Pj.ip to Q\a-\,ii.o is again 

which is precisely the (0, j) entry of the product 

Case 2. j — is similar. 

Case 3. k,j > 1, then the induction hypothesis implies that the sum of paths 
from Pj^i^ to Q\a\ai,k is 

which is precisely the (fc, j) entry of J^-. □ 

When all Bn,k are zero, the level steps of type l^ „ dissapear, and our discussion 
is somewhat related to parts of 

Formula H2.8I) can be also used to calculate the Jacobi coefficients from the 
moments in a relatively simple way. Let a G — {0}. It is convenient to introduce 
the notation i{n) in order to denote the nth letter of the word a (from left to right). 
If \a\ — 2n, then there exists a unique path Po- with corresponding weight 



If |(t| = 2n + l, then there exists a unique path, still denoted Po-, with corresponding 
weight 

In any case, let A4* = A4cr — {Po-}- Also, we introduce the notation: Ai = Ai, and 
for n > 2, 

A - A A®'^ A<BN"-^ 
— ^n^n~l ■ ■ ■ ^1 

Corollary 2.4. The following formulae hold: for n > 1, 



An An — 



Bo,k = Sk, k = l,...,N, 



kl 



— \t\—71/ 



A-^ ■ 



and for n > I, k ~ 1, . . . , N, 
I 



V 



\T\ = \cr\+2=n+l/ 

Due to the fact that An is an upper triangular matrix with strictly positive 
elements on the diagonal, the first relation of the previous result uniquely determine 
An by Cholesky factorization. 

3. Free products 
The set Vn can be viewed as the free product of N copies of Vi: 

= 7^1 ★. . .★Pi = C © (e„>i (B^,^^,,...,^„_,^^„ ^"1 ® • ■ • ® ^l) , 

N times 

where is the set of polynomials in the variable Xi, i = 1, . . . , N, without constant 
term. This remark suggests that the simplest examples of families of orthogonal 
polynomials can be obtained by using free products. Some examples already ap- 
peared in Here we describe a genreal construction. This allows to introduce 
multivariable analogues of all classical orthogonal polynomials. 

The simplest attempt to construct families of orthogonal polynomials on Vn 
would be to consider orthogonal polynomials associated with free products of pos- 
itive functionals. Let (jji and 02 be two strictly positive functional on Vni , respec- 
tively Pn2- Their free product (j> — (j>i -k (j)2 on Vni '*^'Pn2 is defined by 0(1) = 1 and 

0(Pii ...PiJ^ (f>i^{Pi^) . . . (f>i,SPiJ for > 1, «1 «2, ■ • ■, in~l ^ «n, Pit ^ Pn, ' 

and ik G {1,2} for k — 1, . . . ,n. By results in [S], 0], is a positive functional. 
However, as it turns out, is not strictly positive. Thus, consider the case of 
two strictly positive functionals 0i and 02 on Pi. Let K be the restriction of the 
moment kernel Ka, 



K = 



^^02 to the set {1, 12}. Its matrix is then 

01 0i(X2)02(X2) 



02 (X2)01 {Xi) 02 (^2)01 (X/)02 {X2) 



1 

02(^2) 



4>i{xl) 01 (xf) 
01 (x?) 01 (xf) 



1 

02(^2) 



The matrix 



1 1 
1 1 



has rank one, so K is never invertible. 



We can consider another simple way related to free products in order to build 
orthogonal polynomials in several noncommutative variables. Thus, let {^Pn,k\, 
n > 0, k G {l,...,Af}, be N families of orthonormal polynomials on the real line, 
determined by the recursion formulae: 



(3.1) 



X(pn,k{x) = an+l^kfn+l,k{x) + bn,k(Pn,k{x) + {x). 



We introduce polynomials in N noncommutative variables as follows. Any word 
(7 G — {0} can be uniquely represented in the form a = i'l^ . . ■ip'', ii, . . ., 
ip G {1, . . . , N}^ ki, . . kp > 0, and i/ ^ for Z = 1. Then define 

(3.2) ^,{Xl,...,XN)^ipk,,^AX^,)■■■ipkM^^p)■ 

Theorem 3.1. There exists an admissible family A of matrices such that for k = 
1, . . . , and n > 0, 

(3.3) Afc$„ = ^n+lAn+l,k + ^nBn,k + ^n-l^^^fc, 

where = [Vo-]|CT|=n /"'^ n- ^ 0, $_i = 0, and ip^r CLfS- given by H3.2(l . 
Proof. First we prove the result for N — 2. From (|3.1|) we deduce 

bo,i 



ai,i 




so that 
with 

Similarly, 
with 

We see that 



and Bo,i = [ &o,i ] 

X2<^0 = ^'1^1,2 + $0^0,2 

and So, 2 = [ ^0,2 ] 



ai,i 




^1,2 = 



ai.2 




Ai = [ ^1,1 Ai,2 1 = 



ai,i 






01,2 



is upper triangular (actually diagonal) and has the elements on the diagonal > 
(since all a„,fc > 0). The case n > 1 can be delt with in a similar manner. The first 
2"~i words of length n start with letter 1 and have the structure 

where t is a word of lenght k starting with letter 2 (unless it is 0). For k — there 
is exactly one such word, 1", while for < fc < n — 1, there are 2*^"^ such kind of 
words. Using H3.1|l we have that 

Xi(fi^-kT. = Xllplr^-k(pr = XiLpn-kS^Pr 

= ^n-k+l^l^TOn-k+l.l + Vn-kA^rbn-k,! + ¥'n-fe-l,lVTan-fc,l 
= ipi-a-k+l^an-k+l^l + fi^^-k^bn-k,! + fir^-k-irOn-k,!- 



The last 2""^ words of lenght n start with letter 2 and therefore we have for such 
a word a that 

Xupa- = Xnpo^HPa- 

Putting together the above two formulae we deduce that 
where 



O-n-l,! 



^02 



and 



bn,l 



bn-1,1 



1,02 
"n-2,1 



the unspecified entries are all zero. Similarly, we deduce 

for some suitable matrices j4„^2 and Bn,2- Actually, the same proof works for A'' > 2 
and we record here the form of the matrices A„_fe, S„_fe for an arbitrary N. Let 

= {cr G I |(t| = n and a = kr for some r} 

and denote by ttJ! the bijeetion from the set {1, 2, . . . , TV"^-'^} onto WJ! defined 
simply by =the ^th word in Wj^, with respect to the lexicographic order. 

Also, if cr e W" then it has a unique representation a = hPr with r a word that 
does not start with the letter k. Define nk{cr) = p. Now Bn,k is an x matrix 
such that, for l,m £ {1,2, . . . , A'"}, 



(3.4) 



{Bn,k)r 







m 



I 7^ m 



and An^k is an A'" x A'" ^ matrix such that for I e {1,2,. 
{1,2,..'., AT"}, 



A'"-^} and m G 



(3.5) 



{An,k)mi — 



Onfe(7r^(m))-l,fe I = m 

l^m 



Therefore An = [ • • • ^n,jv ] is a diagonal matrix with strictly positive 

diagonal elements, so that A = {An,k,Bm,k | n > 0, m > 0, A; = 1,...,N} is an 
admissible family of matrices. □ 
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